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Abstract: The crosstalk and afterpulsing in Hamamatsu silicon photomultipliers, called Multi- 
Pixel Photon Counters (MPPCs), have been studied in depth. Several components of the correlated 
noise have been identified according to their different possible causes and their effects on the signal. 
In particular, we have distinguished between prompt and delayed crosstalk as well as between trap- 
assisted and hole-induced afterpuising. 

The prompt crosstalk has been characterized through the pulse amplitude spectrum measured at 
dark conditions. The newest MPPC series, which incorporate isolating trenches between pixels, 
exhibit a very low prompt crosstalk, but a small component remains likely due to secondary photons 
reflected on the top surface of the device and photon-generated minority carriers diffusing in the 
silicon substrate. 

We present a meticulous procedure to characterize the afterpuising and delayed crosstalk through 
the amplitude and delay time distributions of secondary pulses. Our results indicate that both noise 
components are due to minority caniers diffusing in the substrate and that this effect is drastically 
reduced in the new MPPC series as a consequence of an increase of one order of magnitude in the 
doping density of the substrate. 

Finally, we have developed a Monte Carlo simulation to study the different components of the 
afterpuising and crosstalk. The simulation results support our interpretation of the experimental 
data. They also demonstrate that trenches longer than those employed in the Hamamatsu MPPCs 
would reduce the crosstalk to a much greater extent. 

Keywords: Avalanche-induced secondary effects; Photon detectors for UV, visible and IR 
photons (solid-state) (PIN diodes, APDs, Si-PMTs, G-APDs, CCDs, EBCCDs, EMCCDs, etc); 
Detector modelling and simulations II (electric fields, charge transport, multiplication and 
induction, pulse formation, electron emission, etc). 
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1. Introduction 

Recent progresses in silicon photomultipliers (SiPMs) have allowed a notable reduction of the cor¬ 
related noise: crosstalk and afterpulsing. However these effects are still major drawbacks of SiPMs. 
Their characterization and modeling are necessary for an appropriate processing and interpretation 
of data as well as for the development of detectors with a further noise reduction. 

Afterpulses originate from parasitic avalanches in a pixel that was previously fired, e.g., by 
an incident photon. They are commonly attributed to the delayed release of carriers that were 
trapped in deep levels in the depletion layer during the previous avalanche []IJ]. There may be 
various afteipulse components due to traps with different lifetimes. Both the trigger probability 
and the charge multiplication of these secondary avalanches are lowered while the pixel voltage is 
recovering from the last avalanche breakdown. Afterpulses distort the single photon counting and 
increase the integrated charge registered for an input light pulse in a time-dependent way. 

The crosstalk is the correlated generation of avalanches in pixels different to the one where 
a primary avalanche took place. Secondary photons emitted in the primary avalanche (see |@] 
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Figure 1. Physical processes that cause correlated noise. Primary avalanches are represented in dark red 
and secondary avalanches in bright red. Optical crosstalk is due to secondary photons reaching directly a 
neighboring pixel (la) or after a reflection on either the top surface (lb) or the bottom surface (Ic) of the 
device. A minority carrier generated in the substrate may diffuse to the primary pixel (2a) or to a neighboring 
one (2b) and produce an afterpulse or delayed crosstalk, respectively. Trap-assisted afterpulsing (see text) is 
not illustrated in the figure. The electric held intensity as a function of depth is shown on the right. 

and references therein) can reach neighboring pixels, either directly or after a reflection, and trig¬ 
ger avalanches there almost simultaneously to the primary one (figure |T]). Hence, this effect is 
called optical crosstalk. Secondary avalanches also emit photons that can trigger more simultane¬ 
ous avalanches in other pixels, giving rise to high-amplitude pulses even at dark conditions. The 
crosstalk degrades the excellent photon counting resolution of SiPMs, unless the crosstalk multi¬ 
plicity is known and properly accounted for [Q]. 

Besides, secondary photons absorbed in non-depleted regions of the device, especially the 
silicon substrate, produce carriers that can diffuse to the depleted region of a pixel and trigger a 
secondary avalanche [Q]. In the case of the p on n structure shown in figure [I], only holes originated 
from the n-type substrate are accelerated upwards by the electrical field in the junction and able 
to trigger an avalanche. The delay time of this secondary avalanche is dominated by the canier 
diffusion time, which relies on the depth at which the secondary photon was absorbed and the 
doping density of the substrate. Depending on whether the carrier reaches the primary pixel or any 
other one, the event is classified as afterpulse or crosstalk, respectively. These two components 
are undistinguishable at long delay time, but afteipulses with delay times shorter than the pixel 
recovery time have lower amplitude. On the other hand, since crosstalk avalanches are not affected 
by the primary pixel recovery, they may have arbitrarily short delay times. If the delay time is even 
shorter than the time resolution of the detection system, the output signal will be basically the same 
as if the secondary photon had been absorbed within the depletion region of the pixel. According 
to that, crosstalk events are classified as prompt or delayed. 

In this work, we focused on the SiPMs fabricated by Hamamatsu Photonics K.K. [0. This 
company named its product Multi-Pixel Photon Counter (MPPC), which is a trademark. Hama¬ 
matsu has released several series of MPPCs based on the p on n “reverse junction” structure of 
figure [I], which was designed to optimize the photon detection efficiency around 400 nm |0. The 
S10362-11 and S10362-33 series (not marketed now) are MPPCs with active areas of 1 mm x 1 mm 
and 3 mm x 3 mm, respectively. They were fabricated with pixel pitches of 25, 50 or 100 /im. The 
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S12571 and S12572 series were produced with the same active areas as the former S10362-11 and 
S10362-33 series, but in addition to the 25, 50 and 100 /im pixel pitch types, MPPCs with pixel 
pitches of 10 and 15 /^m were produced, offering larger dynamic range and shorter recovery time. 
A high fill factor (i.e., the ratio of the active area of a pixel to its total area) ranging from 33% 
for the 10 /im pixel pitch type to 78% for the 100 /im pixel pitch type was achieved. The most 
prominent progress with respect to the older series is a drastic reduction of afterpulsing, which 
was attained thanks to the use of improved materials and wafer process technologies according 
to Hamamatsu’s reports. Recently, it has been released the SI3360 series incorporating optical 
barriers, called trenches, between pixels that drastically reduce crosstalk too. Low crosstalk and 
afterpulse MPPCs with active areas of 1.3 mm x 1.3 mm, 3 mm x 3 mm and 6 mm x 6 mm, and 
with pixel pitches of 25 and 50 /xm are available. We characterized experimentally and modeled the 
afterpulsing and crosstalk for these three series of MPPCs, studying the different noise components 
according to their physical causes. Solutions for a further reduction of these effects are proposed. 

The paper is structured as follows. The experimental method is described in section The 
results on prompt crosstalk are given and compared with model predictions in section ^ The after¬ 
pulsing and delayed crosstalk are characterized and modeled in section 0. The possible causes of 
correlated noise and their solutions are analyzed by means of Monte Carlo simulations in section^ 
The conclusions are in section 0. 

2. Experimental method 

In this work, all measurements were carried out at dark conditions and room temperature; the 
detector temperature ranging from 25 to 30°C. Pulses produced by thermal effect, called dark 
counts, were registered at a low rate of < 1 MHz. This allowed a characterization of the correlated 
noise on a single primary avalanche basis and almost free from contamination due to coincident 
uncorrelated dark counts. Our technique is based on a waveform analysis to identify individual 
pulses through their leading edges. We focused on the MPPCs with the smallest active area (i.e., 
1 mm X 1 mm and 1.3 mm x 1.3 mm), because they produce narrower and higher pulses, although 
some 3 mm x 3 mm active area MPPCs were also characterized. 

We employed a Hamamatsu C12332 driver circuit for MPPCs 0. It included a regulable 
power supply, a high-speed operational amplifier and a temperature sensor. Output signals were 
registered by a digital oscilloscope (Tektronix TDS5032B) and stored in binary files for later analy¬ 
sis. The oscilloscope was programmed to trigger randomly 1000 times per run with a time window 
of 3.2 ms and a sampling rate of 2.5 GS/s. A run contained typically around 10^ pulses. Tempera¬ 
ture variations during a run were less than 0.5°C and were neglected in the analysis, but run-to-run 
temperature variations were corrected using the temperature coefficient of the breakdown voltage 
provided by Hamamatsu. 

The data analysis procedure was described in detail in [@] and only a few upgrades were made 
since then. A C-i-i- code was developed to read the files generated by the oscilloscope’s software 
and perform a waveform analysis. It deconvolves the signal assuming an exponential response 
function with a time constant depending on the MPPC type. The detectors of the S13360 series 
were found to produce pulses with a spike component of a few ns followed by a slow component 
of tens of ns (see, e.g., 0). The deconvolution only removes the slow component in this case. A 
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pulse is identified when the deconvolved signal surpasses a discrimination threshold that is tuned 
for each MPPC and bias voltage. In this way, pulses are separated with a resolution better than 
10 ns. Furthermore, for every identified pulse, an algorithm searches for double-peak structures in 
the signal over the threshold revealing very close pulses that cannot be separated. The resolution 
to identify these piled-up pulses is estimated to be 3 ns, which determines the minimum delay time 
to distinguish delayed from prompt crosstalk. 

The pulse amplitude and time are measured at the maximum of the deconvolved pulse. In 
addition, the code calculates the pulse height in the original signal after subtracting the baseline 
by means of an exponential extrapolation of the signal in a small region prior to the pulse edge. 
The original pulse height has a lower statistical uncertainty because it is less affected by electronic 
noise than the deconvolved pulse maximum, but the former cannot be applied to a pulse that is 
closer than 20 ns to the preceding one. 

We studied the dependence of the correlated noise on the MPPC gain, which is proportional 
to the overvoltage, i.e., the difference between the bias voltage and the breakdown voltage. To do 
that, the breakdown voltage of each MPPC was determined by zero extrapolation of the separation 
between the first and second peaks in the pulse amplitude spectrum versus bias voltage [@]. Both 
the deconvolved pulse amplitude and the original pulse height gave the same breakdown voltage 
within uncertainties. The proportionality constant that relates overvoltage and gain was provided 
by Hamamatsu. 

3. Characterization of prompt crosstalk 
3.1 Pulse amplitude spectra 

We characterized the prompt crosstalk from either the deconvolved pulse amplitude spectrum or 
the pulse height spectrum, although the latter had a better resolution in most cases. Only pulses 
that were separated from the previous pulse by more than 120 ns were selected to calculate the 
spectrum. This prevented including pulses originated from pixels that were still recovering, as 
they have lower gain and therefore lower probability to induce crosstalk. Cuts were also applied 
to remove piled-up pulses and false pulses due to electronic noise exceeding the discrimination 
threshold. 

The relative areas of the peaks in the amplitude spectrum give the probability distribution 
of the prompt crosstalk multiplicity, where the k-th peak corresponds to a multiplicity k-1. The 
overall crosstalk probability e is thus defined by one minus the relative area of the first peak. More 
precisely, a small correction was made to account for coincidences of uncorrelated events in a time 
interval of 3 ns (< 0.001 probability), as was described in |^. 

The spectra obtained for the S10362-11 and S12571 series show clear Gaussian peaks (fig¬ 
ure 0). Thus, a multi-Gaussian tit was performed to correct for peak overlapping for those MPPCs. 
On the other hand, for the new S13360 series, with trenches, the spectra exhibit asymmetrical peaks 
for crosstalk multiplicity > 1, especially when using the deconvolved pulse amplitude (figure 
This is due to pulses formed by two or more avalanches adding up in a time interval less than 3 ns, 
since these events have to pass the cut of piled-up pulses, but asynchronously enough to give a 
pulse amplitude lower than the expected one for perfectly synchronous avalanches. The above- 
mentioned spike component of the pulses enhances this effect, which is much more noticeable in 
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Figure 2. Pulse height spectrum obtained for the S12571-050C MPPC at dark conditions. Both a high 
probability and multiplicity of prompt crosstalk are observed. Our waveform analysis allowed us to resolve 
at least 16 Gaussian-like peaks. 



Pulse height (mV) Deconvolved pulse amplitude (mV) 

Figure 3. Pulse height (left) and deconvolved pulse amplitude (right) spectra obtained for the S13360- 
1350CS MPPC at dark conditions. The probability of prompt crosstalk is much smaller than in the S12571- 
050C MPPC despite the higher gain. The peaks are asymmetrical because some crosstalk avalanches are not 
perfectly synchronous with the primary one. This effect is much more noticeable in the deconvolved signal 
because of the shortening of pulses. 


the deconvolved signal because of the shortening of pulses. The trenches clearly reduce the prompt 
crosstalk, but this lack of synchronism suggests that a significant fraction of the prompt crosstalk 
remaining is induced by minority carriers generated at very low depth in the substrate, that is, close 
to the depletion region of neighboring pixels, as will be demonstrated in section ^ by Monte Carlo 
simulation. Coincidences of uncorrelated pulses cannot justify it, because the dark count rate is 
extremely low. 

3.2 Overall prompt crosstalk probability 

The overall prompt crosstalk probability e versus gain is shown in figure 0 for mosf of fhe defecfors 
fhaf were characferized. The uncerfainfies, which are basically sfafisfical (see [Q] for defails), are 
less or around 1% af moderafe and high gain, excepf for fhe S10362-11-025C, S13360-1325CS 
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Figure 4. Overall prompt crosstalk probability as a function of gain for various MPPCs. Symbols represent 
experimental data and lines represent the fits of eq. (|3.1D to data. 


and S13360-3025CS MPPCs, for which the uncertainties are around 10%. It can be seen that £ 
grows more than linearly with gain. The number of emitted secondary photons is proportional to 
the avalanche carrier multiplication, but the probability of a photon to trigger an avalanche in a 
neighboring pixel also increases with overvoltage. On the other hand, £ must saturate at high gain. 
According to that, the following function was fitted to data: 


£ (X) = 1 — exp 



(3.1) 


where X can represent either the gain or the overvoltage, and X is a MPPC constant in the cor¬ 
responding units. The index a is positive and accounts for the increase of the trigger probability 
with overvoltage assuming a power law. In principle, the overvoltage dependence of the trigger 
probability can be extracted from the curve of the photon detection efficiency (PDE) versus over- 
volfage, which is provided by Hamamafsu. This curve can be described by a power law wifh 
exponenf around 0.5, buf wifh a positive overvolfage offsef fhaf increases wifh decreasing pixel 
size. Nonefheless, fhe above funclion fifs well fo dafa, as shown in fhe figure. 

The filled K and a values, as well as fhe proporlionalify consfanl belween gain G and overvoll- 
age AV, are lisled in fable |T]. The K and a paramefers are correlaled wifh each olher and Iherefore 
Iheir uncerfainlies may be large. Neverlheless, fhe filled £ values agree wifh fhe measured ones 
wilhin experimenfal uncerfainlies. 

Our resulls are summarized in fhe following poinls: 


1. The prompf crosslalk increases rapidly as fhe pixel size decreases keeping fhe same gain. 
This is because secondary photons have more chances fo escape a small pixel. 

2. The a index also increases wifh decreasing pixel size. The cause is nol clear, buf Ibis vari- 
alion is consisfenf wifh fhe above-menfioned shiff of fhe PDE curve towards positive over- 
volfages. 
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Table 1. Proportionality constant between gain G and overvoltage AV, and the fitted K and a values to 
prompt crosstalk data for all the MPPCs that were characterized. 


MPPC 

Hi^(xlOVV) 

^(V) 

a 

S10362-11-025C 

0.133 

10.57 

1.20 

S10362-11-050C 

0.580 

5.40 

0.65 

S10362-11-100C 

2.360 

2.79 

0.45 

S10362-33-100C 

2.360 

2.56 

0.62 

S12571-010C 

0.030 

12.81 

1.71 

S12571-025C 

0.158 

7.37 

0.71 

S12571-050C 

0.506 

4.98 

0.60 

S12572-050C 

0.506 

4.03 

0.56 

S12571-100C 

2.068 

2.30 

0.52 

S13360-1325CS 

0.129 

103.24 

0.82 

S13360-3025CS 

0.129 

32.23 

1.50 

S13360-1350CS 

0.571 

49.38 

0.58 

S13360-3050CS 

0.571 

39.42 

0.40 


3. The S12571 series have a prompt crosstalk probability higher than the old S10362-11 series, 
especially for the 25 /tm pixel pitch type. In fact, Hamamatsu reports a fill factor of 65% for 
the S12571-025C MPPC in contrast to 30.8% for the S10362-11-025C MPPC, although the 
fill factors of the 50/im and 100 /im pixel pitch types are very similar in both series. 

4. The prompt crosstalk probability for the S13360-1325CS and S13360-1350CS MPPCs is 
lower by a factor of around 80 and 35, respectively, than for the S12571-025C and S12571- 
05OC ones. 

5. The 3 mm x 3 mm active area MPPCs have somewhat higher prompt crosstalk probabilities 
than their small active area counterparts, which have a lower number of pixels. 

3.3 Prompt crosstalk multiplicity 

The probability distribution of the prompt crosstalk multiplicity for a given e value depends on the 
number of pixels that can be reached by secondary photons from the primary pixel and on crosstalk 
cascading, i.e., secondary avalanches inducing further avalanches. This was modeled in a previous 
work 10], where exact analytical formulae were derived from simple geometrical assumptions. 
For the MPPCs that were characterized in the present work, two models were found to fit data 
in general, namely that either the 4 nearest neighbors or the 8 nearest neighbors of each fired 
pixel have equal probability of prompt crosstalk, whereas the probability is zero in the remaining 
pixels. In figure ^ we show the relative difference of the multiplicity distributions predicted by both 
models with respect to the experimental one for 50 /im pixel pitch MPPCs and up to multiplicity 4. 
The model predictions were normalized so that the probability of multiplicity 0 equals 1 — e. Also 
shown in the figures are the results for two other analytical models proposed in [0 , which represent 
limit situations regarding the number of neighboring pixels subject to prompt crosstalk (see |@] 
for details). An empirical formula proposed in HTHH proved to describe the crosstalk multiplicity 
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Figure 5. Relative differences of the prompt crosstalk multiplicity distributions predicted by different an¬ 
alytical models with respect to the experimental one for 50 fim pixel pitch MPPCs. Results suggest that 
prompt crosstalk affects at least the 8 nearest neighboring pixels, but with a higher probability for the 4 near¬ 
est ones. Uncertainties are large for the S13360-1350CS MPPC because the prompt crosstalk probability is 
very low. 


distribution for the S10362-33-050C MPPC, but we did not included it in the comparison because 
it uses a fitting shape parameter and has no clear geometrical interpretation. 

The following features were observed in the comparison of data with the model predictions: 


1. For a given MPPC series, the model that agrees with data the best is independent of the bias 
voltage and the size of the active area (see also |0). This indicates that the prompt crosstalk 
is mostly determined by the pixel structure. 

2. The larger the pixel size, the lower the number of neighboring pixels subject to prompt 
crosstalk according to our models. Although this trend is subtle, it is the expected one for an 
exponential attenuation of secondary photons. Small differences were also observed between 
the S10362-11 and S12571 series for the same pixel pitch. 

3. The agreement of models with data for the S13360 series is compatible with the other series. 
However, the uncertainties are large due to the low probability of prompt crosstalk and no 
conclusions can be drawn. 
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Figure 6. 2D histogram of deconvolved pulse amplitude versus elapsed time from the last pulse obtained 
for the S10362-11-050C MPPC at an overvoltage of 2.1 V (G = 1.2 • 10®). Data are distributed in several 
horizontal bands corresponding, from bottom to top, to prompt crosstalk multiplicity 0, 1, 2, ... The low- 
amplitude branch of each band is made up by afterpulses. The grey square indicates the secondary pulses 
selected to obtain the distribution ^noAp(0- 

4. Characterization of afterpulsing and delayed crosstalk 
4.1 Procedure 

In a previous work |p, we presented a procedure to characterize the afterpulsing and delayed 
crosstalk through their delay time and amplitude distributions. The procedure has been upgraded 
and it is described in detail here. 

The first step was to obtain the 2D histogram of the deconvolved pulse amplitude versus the 
elapsed time from the last pulse for each MPPC and bias voltage. For the sake of clarity, a pulse 
contributing to this histogram will be referred to hereafter as a secondary pulse, and its preceding 
pulse giving the origin of time, as the primary one, although they are not necessarily correlated 
with each other. To obtain the histogram, we only selected secondary pulses of which the primary 
pulse has an amplitude within the range of values corresponding to a single avalanche (i.e., no 
prompt crosstalk) and passes the cut of piled-up pulses. The primary pulse was also imposed to 
be separated by at least 500 ns from previous ones to prevent contamination of secondary pulses 
correlated with them. An example of such a histogram for the S10362-11-050C MPPC is shown 
as a color map in figure ^ The undermosf horizonfal band is made up by secondary pulses wifhouf 
prompf crossfalk, and the upper bands correspond to prompt crosstalk multiplicities equal or greater 
than 1. The low-amplitude branch of each band is due to afterpulses that are produced while the 
pixel is still recovering from the primary pulse. Occurrence of two consecutive uncorrelated dark 
events in the same pixel also contributes to this branch, but it is very unlikely. Full-amplitude pulses 
at short time are delayed-crosstalk and uncorrelated dark events (in different pixels). 

In a second step, the amplitude-time correlation of afterpulses was utilized to characterize 
the recovery function of the MPPC. The 2D histogram was integrated in intervals of 2 ns limited 
to the first band of secondary pulses, obtaining a collection of amplitude histograms, each one 
at a different stage of the pixel recovery. This is shown in figure [3, where fwo Gaussian-like 
componenfs, corresponding fo afferpulses and full-amplifude pulses, are disfinguished for every 
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Figure 7. Pulse amplitude distributions extracted from the 2D histogram of hgure ^by integrating in 2 ns in¬ 
tervals. The distributions are restricted to single-avalanche pulses (no prompt crosstalk). A double-Gaussian 
fit was performed to characterize the pixel recovery and the losses of afterpulses with amplitude below the 
discrimination threshold (Athr = 14 mV in this case). 


time interval. We observed that the variance of both Gaussian components increases linearly with 
increasing gain, as expected from the Poissonian nature of the avalanche carrier multiplication. 
Thus, a double-Gaussian fit was performed with the following restrictions: a) the mean of the 
component of full-amplitude pulses was set equal to the mean value Ai of the first peak in the 
deconvolved pulse amplitude spectrum described in section b) the variance of both Gaussian 
components was forced to be proportional to the mean, where the proportionality constant C was 
taken as a common fitting parameter for all the time intervals, and c) the mean of the component of 
afterpulses was assumed to vary with the time interval as 


=Ai • 6{tk — to) ■ 


1 — exp 




(4.1) 


where is the central value of the ^-th interval, 6{t) is the Heaviside step function, tj-ec is the 
recovery time and to is a positive time offset that was needed to fit data for some MPPCs. This 
offset could arise from the fact that the pixel voltage may drop below the breakdown voltage m- 
As expected, the fitted tree, C and to values showed no dependence on bias voltage. Average free and 
fo values for every MPPC are shown in table 

In a third step, the 2D histogram was integrated over amplitude and normalized to obtain the 
probability density distribution of delay time P{t), including all types of the secondary pulses. 
Although some pulses with very short delay time were able to be resolved, we set a minimum 
delay time fmin that defined fhe effeefive time resolution (tmin was eifher 6 ns or 10 ns, depending 
on fhe MPPC). The resulting P(f) disfribufion is represenfed by black open circles in figures and ^ 
for fhe S10362-11-050C and S13360-1350CS MPPCs, respectively. In addition, we obfained fhe 
probabilify densify disfribufion PnoAp(0 where fhe afferpulses were removed for delay time shorfer 
fhan a cerfain value fmax depending on fhe MPPC (red open circles in figures ^ and To do 
fhaf, we firsf obfained a disfribufion P°oap( 0 selecting only full-amplilude pulses of fhe undermosf 
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Figure 8. Delay time distribution of secondary pulses for the S10362-11-050C MPPC extracted from figure 
In the left-hand plot, two components of trap-assisted afterpulses with different lifetimes were fitted to 
data. The fit is good, but the high density of full-amplitude pulses at short time (red open circles) cannot be 
justified by uncorrelated dark counts only. In the right-hand plot, only the components of afterpulsing and 
delayed crosstalk due to minority carriers diffusing in the substrate were fitted to data. Both the P{t) and 
PnoAp(0 distributions are properly described by this model. 
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Figure 9. The same as figure || for the S13360-1350CS MPPC. The correlated noise decays much faster 
than for the S10362-11-050C MPPC. 


band of the 2D histogram, that is, delayed-crosstalk and uncorrelated dark events without prompt 
crosstalk, as illustrated by the grey square in figure Then, the PnoAp(0 distribution was estimated 
by 


PnoAp(0 = (l+£)-<Ap(0, (4.2) 

where the factor 1 + £ adds the probability of secondary pulses with prompt crosstalk, but still 
excluding afterpulses with or without prompt crosstalk. 

The last step was to find an appropriafe combination of noise components that describes the 
experimental distributions P{t) and PnoAp(0. The P{t) distribution can be expressed as 


P{t) = exp 


1 

>1 


^min 



Pdc + X1/K0 


(4.3) 


where Pdc is the dark count rate and /, (t) is the density distribution of secondary pulses of type 
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i directly correlated with the primary pulse, i.e., excluding cascades. To derive the above equa¬ 
tion, we assumed that the number of secondary pulses in a given time interval follows a Poisson 
distribution. Likewise, the PnoAp(0 distribution is expressed as 


^noAp(0 = exp 


-^DC • (t — tmin) — f fcri^) ' 

^min 


[■^DC+/ct( 0 ] ! 


(4.4) 


where /ct (0 is the density distribution of delayed-crosstalk pulses. 

Note that P{t) is proportional to exp(—7 ?dc • 0 at long delay time, because the correlated noise 
extinguishes with time. We used it to determine 7 ?dc by means of an exponential fit for delay time 
greater than 1 /is. As expected, the fitted Rue values were proportional to overvoltage and were 
lower than 1 MHz in most of the cases, with the new MPPC series having the lowest Rue values. 
Even if a /is-life component of correlated noise were present, it would be undistinguishable from 
the uncorrelated dark noise. 

To model the correlated noise components, we proceeded in the following way. Each compo¬ 
nent of trap-assisted afterpulses was parameterized by 


/trap(0 ~ ^trap • exp I | •/I(t) • k(t) , (4.5) 

\ hrap / 

where the coefficient Qrap only depends on overvoltage and Ttrap is the trap lifetime. The factor 
/l(t) accounts for the recovery of the trigger probability in the primary pixel and it is assumed to 
be 


Hit) 


Ai ’ 


(4.6) 


where A(t) is given by (|4.![). The factor k{t) is the fraction of afterpulses that surpass the discrimi¬ 
nation threshold Athr at a given delay time. Note that afterpulses with prompt crosstalk have always 
amplitudes above Athr- Taking into account that the probability of an avalanche to induce prompt 
crosstalk is proportional to its carrier multiplication, the time-varying fractions of afterpulses with 
and without prompt crosstalk can be estimated by £ • /i(t) and 1 — £ • /i(t), respectively. Now, let 
ko{t) be the fraction of afterpulses without prompt crosstalk that surpass the threshold. Assuming 
a Gaussian distribution of afterpulse amplitude at given delay time, as explained above, ko{t) is 


koit) = 


Uthr y^2-K-C-A{t) 


■ exp 


where C is already known. Therefore, k{t) is given by 


{x-A{t)f 

2-C-A{t) 


(4.7) 


k{t) = (l-£-/i(t))-ko(f)-h£-/l(t). (4.8) 

We also considered the components of afterpulsing and delayed crosstalk induced by minority 
carriers diffusing in the substrate. As will be explained in section ^ a Monte Carlo simulation 
was performed to evaluate the diffusion time distribution for minority carriers reaching either the 
primary pixel or any other pixel. In both cases, the distribution was found to be well described 
by the product of a power law and an exponential decay function related to the minority carrier 
recombination. Therefore, we obtained 
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/ap(0 = Cap • C • exp 



(4.9) 


and 


/ct( 0 =CcT-?^-exp 


"Z^bulk 


(4.10) 


where the a and b exponents are both negative and Tbuik is the minority carrier lifetime in the silicon 
substrate, also called bulk lifetime. Notice that the delayed crosstalk is not affected by the /i(t) and 
k{t) factors. 

The number of trap-assisted afterpulse components, if any, is not known a priori. Besides, 
many different combinations of noise components may fit to the experimental time distributions. 
Therefore, we opted for making two simple hypotheses. In the first one (HI), we considered 
exclusively two trap-assisted afterpulse components: “slow” and “fast” (see, e.g., EUl); four fitting 
parameters were defined: Csiow; fsiowj Cfast and Tfast; and fhe PnoAp(0 disfribufion was nof used in 
fhe fif. In fhe second one (H2), only afferpulsing and delayed crossfalk relafed fo caiTier diffusion 
were considered; fhe fiffing parameters were Cap, Cct and Tbuik; we sef a = — 1 and b = —0.5, 
according fo our simulafion resulfs (see secfion fo reduce fhe number of degrees of freedom; 
and we puf fhe exfra resfricfion fhaf fhe sum of fhe dark-counf and delayed-crossfalk componenfs 
has fo fif fo /’noAp(0- The resulfs for bofh hypofheses are discussed in fhe nexf secfion. 


4.2 Results 


The hypofhesis HI was found fo describe successfully fhe shape of P{t) for all fhe MPPCs and 
bias volfages. As an examples, fhe besf fifs fo P{t) for fhe S10362-11-050C and S13360-1350CS 
MPPCs are shown in fhe leff-hand plofs of figures ^ and ^ respecfively, where fhe sfacked areas 
represenf fhe fhree componenfs considered: dark counfs and bofh fypes of frap-assisfed afterpulses. 
However, HI has two deficiencies fhaf make if incompafible wifh dafa. In fhe firsf place, a signifi- 
canf fraction of pulses fhaf populate P{t) af shorf delay fime are full-amplifude pulses fhaf belong 
fo PnoAp(0 and cannof be justified solely by dark counfs. If is clear fhaf fhe delayed crossfalk is 
non-negligible. In fhe second place, fhe tiffed Tsiow and Tfast values, which represenf supposedly 
fhe lifefimes of fwo fypes of fraps in fhe deplefed region, vary nofably from MPPC fo MPPC even 
wifhin fhe same series. For insfance, we obfained Tgiow = 70 ns and Tfast = 10 ns for fhe S10362- 
11-050C MPPC, buf Tsiow = V ns and Tfast = 4 ns for fhe S13360-1350CS MPPC. 

For fhe hypofhesis H2, fhe fif fo P{t) was as good as for HI, despite fhe facf fhaf fhe number 
of fiffing parameters was reduced from 4 fo 3. In addition, a good fif fo PnoAp(0 was achieved by 
including fhe componenf of delayed crossfalk. The besf fifs for fhe S10362-11-050C and S13360- 
1350CS MPPCs are shown in fhe righf-hand plofs of figures and 0, respecfively. Note fhaf fhe 
power-law factors in equations (|4.9|) and (|4.10|) properly accounf for fhe non-exponenfial fails of 
P{t) and PnoAp(0- The small inaccuracies of fhe fifs (< 10%) are likely due fo fhe assumptions 
made on fhe indexes a and b as well as on fhe jU(t) and k{t) functions, which defermine fhe fasf 
rising and fhe peak of P{t). Correlafions were also found befween all fhese parameters and fhe 
fiffing paramefers Cap, Cct and Tbuik- Adding a frap-assisfed afterpulse componenf would improve 
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Table 2. Fitting parameters of the model for afterpulsing and delayed crosstalk. The integrals defined in 
eq. (l4.12| ) calculated for G = 10^ are given in the last two columns. 


MPPC 

free (ns) 

fo (ns) 

Tbuik (ns) 

^AP (V) 


S (ns-"'^) 

PkP 

VcT 

S10362-11-025C 

5.0 

0.6 

73.0 

11.78 

1.40 

0.060^ 

0.74 

0.25 

S10362-11-050C 

12.5 

0.0 

120.0 

4.86 

1.05 

0.021 

0.25 

0.04 

S10362-11-100C 

36.0 

4.3 

244.0 

2.26 

0.67 

0.015 

0.09 

0.02 

S10362-33-100C 

43.7 

2.4 

395.0 

1.88 

1.47 

0.017 

0.04 

0.01 

S12571-010C 

4.4 

0.0 

6.4 

25.96 

1.71 

0.140'’ 

0.88 

0.43 

S12571-025C 

8.5 

0.0 

10.8 

11.33 

1.15 

0.060 

0.15 

0.05 

S12571-050C 

20.5 

1.0 

12.3 

6.24 

0.89 

0.021 

0.04 

0.01 

S12572-050C 

20.4 

0.0 

9.6 

5.36 

1.00‘= 

O.O 2 U 

0.04 

0.01 

S12571-100C 

56.1 

1.1 

11.0 

2.04 

0.55 

0.006 

0.01 

<0.01 

S13360-1325CS 

16.7 

0.0 

9.4 

15.90 

1.07 

0.044 

0.07 

0.02 

S13360-3025CS 

12.4 

0.0 

10.8 

21.64 

1.00‘= 

0.076 

0.06 

0.03 

S13360-1350CS 

29.0 

1.4 

8.5 

6.09 

1.13 

0.022 

0.01 

<0.01 

S13360-3050CS 

28.5 

0.0 

11.4 

9.90 

1.04 

0.030 

0.01 

<0.01 


^ Value taken from the S12571-025C MPPC. 
Value obtained by simulation. 

Non-fitted value. 


Value taken from the S12571-050C MPPC. 


the fit, but this does not seem to be necessary to describe the shape of P{t). Some authors, who car¬ 
ried out a characterization of the afterpulsing and delayed crosstalk in other SiPMs, attributed the 
afterpulses to trap-assisted processes only [ ]T31 , [1^ . However, the occurrence of delayed crosstalk 
induced by minority carriers diffusing in the substrate should be along with a component of after¬ 
pulsing. Consequently, we adopted the hypothesis H2 to interpret our data. 

The fitted Tbuik values are a few hundreds of ns for the old S10362 series and around 10 ns for 
the newer ones. This seems to be the main cause of the afterpulse reduction in the new MPPCs, 
since the lower the bulk lifetime of the substrate, the lower the number of minority carriers reaching 
the depleted region of pixels. Using data from []T^, these Tbuik values translate into an increase of 
the substrate doping density from 3-6 -10*^ cm^^ for the S10362 series to 2 • lO'^ cm^^ for the 
newer ones. The recovery time also increases from the old to the new series, contributing to the 
afteipulse reduction too. These results are consistent with the findings of []T^. 

The Cap and Cct coefficients were found to have a power-law dependence with gain, but their 
ratio keeps approximately constant. Since the two noise components are closely related to each 
other, we preferred to characterize them by taking average values of Tbuik and ^ for each MPPC 
(table 0) and refit Cap as the only parameter depending on gain. The fitted Cap values (dimension¬ 
less units) for most of the MPPCs are shown in figure |T^. MPPCs belonging to different series but 
having the same pixel size have very close Cap values at a given gain, despite the large differences 
in their total amount of afterpulses. Similarly to the prompt crosstalk probability (figure Cap 
increases with decreasing pixel size, although the cause of this behavior is not well understood 
yet. The ^ ratio also shows a strong correlation with the pixel size, which was studied by Monte 
Carlo simulation (section For some MPPCs where the components of afterpulsing and delayed 
crosstalk could not be separated, we used either the ^ value measured for other MPPC with same 
pixel size or the value obtained by simulation. 
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Figure 10. Fitted Cap values as a function of gain for the MPPCs with small active area. The smaller the 
pixel size, the greater the Cap value. The lines represent power-law fits to data. 


The dependence of Cap on gain or overvoltage is parameterized in a similar way to e in 
eq. O, with the difference that Cap is not a probability and thus it may take values greater 
than unity: 

/ X \ 

The best fits to data are shown in figure m and fhe ^ap and j3 values are given in fable ^ As 
expecfed, fhe j3 indexes are similar fo fhe a values obfained for e (fable [I]). In some cases where 
fhe dafa were scarce, we sef j3 = 1 . 

To quanfify fhe correlafed noise of a MPPC, fhe mean number of secondary pulses of each fype 
per primary avalanche can be evaluafed by infegrafion of equations (|4.9|) and (|4.1U| ). The following 
infegrals were calculafed for a reference gain of 1 • 10^: 


Nap =Cap ■ [ 7 • exp (-• ju(t) • dl 

JSns t y Tbulk / 

Nct =Cct• / T .exp (-(4. 12 ) 

JSns \/t \ fbulk / 

where we assumed fhaf fhe /ap(0 and fci{t) functions can be exfrapolafed down fo 3 ns. Nofe 
also fhaf fhe k{t) factor in eq. ( |4.9| ) has been removed here fo include all fhe afferpulses, regardless 
fhe discriminafion fhreshold. The resulfs of fhese infegrals are given in fable |I]. Bofh Nap and Act 
decrease by a factor of 5 or more from fhe old S10362-11 series fo fhe S12571 series. A furfher 
reduction of a factor of around 2 is observed in fhe S13360 series, buf if is nof as drastic as fhe 
reduction of prompf crossfalk. Trenches seem fo have no significanf effecf on delayed crossfalk. 

We made a furfher fesf of hypofhesis H2. The analysis was repeated for fhe S10362-11-050C 
and S13360-1350CS MPPCs buf selecting primary pulses wifh prompf crossfalk mulfiplicify 1, 
i.e., fwo simulfaneous primary avalanches. As expecfed, Tbuik mainfained fhe same value and Aap 
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increased approximately by a factor of 2; more precisely 2.08 for the S10362-11-050C MPPC 
and 2.03 for the S13360-1350CS one. On the other hand, the increase in A^ct was only 35% and 
37%, respectively. This could be justified by the fact that crosstalk between the two primary pixels 
during their recovery gives rise to low-amplitude pulses that would be classified as afferpulses. A 
generalization of our model fo be applied fo an arbifrary disfribufion of primary avalanches (e.g., 
friggered by a lighf pulse) should include cascades of correlafed pulses and fhe recovery of every 
fired pixel, nof only fhe primary ones. Implemenfing fhese feafures would also need fo model whaf 
neighboring pixels are subjecf fo delayed crossfalk for each fired pixel, as done for fhe prompf 
crossfalk. 

5. Monte Carlo simulation 

We developed a Monfe Carlo simulation of fhe prompf crossfalk, afferpulsing and delayed crossfalk 
for MPPCs. The goal of fhis simulafion is fwofold. In fhe firsf place, if has provided us wifh 
fhe informafion needed fo parameferize fhe time disfribufions /ap(0 and /ct(0 for modeling fhe 
afferpulsing and crossfalk due fo minorify carriers diffusing in fhe silicon subsfrafe (secfion O)- 
In fhe second place, if allowed us fo sfudy some aspecfs of fhe differenf componenfs of crossfalk 
and afferpulsing, as fheir dependence on geomefrical and physical paramefers and fhe impacf of 
frenches. 

The simulafion includes fhe emission of secondary photons from a primary avalanche, fheir 
absorpfion in fhe differenf parts of fhe device and fhe diffusion of phofon-generafed minorify carri¬ 
ers in fhe subsfrafe. Bofh phofons and carriers can frigger new avalanches, buf fhe phofons emiffed 
by fhese secondary avalanches are nof simulafed, i.e., cascades are ignored. 

5.1 Input physical parameters 

The wavelengfh specfrum of secondary phofons was faken from [@], which is reproduced in fig¬ 
ure |T^ in relafive unifs. For simplicify, fhe disfribufion was sampled in 60 discrefe wavelengfh 
values befween 450 and 1620 nm. This spechum was measured for a MPPC of fhe S10362-11 
series, allhough if is expecled fo be valid for any SiPM. The absorption lenglh of silicon al room 
lemperalure as a function of wavelengfh was faken from [1T511 . 

The diffusion conslanf of minority carriers (holes) in fhe n-fype subsfrafe was varied wilhin a 
large range from 0.1 fo 1 /im^/ns, buf we generally used a reference value of 0.4 /rm^/ns corre¬ 
sponding to a doping density equal or higher lhan 5 • 10^^^ cm^^ [jl^ . 

Some basic geomefrical parameters of fhe MPPCs, as fhe pixel pilch Lpix, fhe fill faclor F and 
fhe size of fhe aclive area are provided by Hamamalsu. Besides, we carried oul measuremenfs on 
some MPPCs of fhe differenf series to verify and complefe fhis informafion. Firsf, we decapsu- 
lated and recovered fhe silicon dice infegraled in fhe defecfor and measured fhe surface value of 
ifs geomefrical paramefers (e.g., die size, pixel size and separalion befween adjacenl pixels) using 
an oplical microscope lhal achieves a magnificalion value of 1500. Then, by means of a grinding 
and lapping mefallographic process, we made several cross-seclions in differenf areas of inferesl 
of fhe analyzed sample, which were observed fhrough bofh and optical microscope and a scanning 
eleclron microscope (SEM). This allowed us fo determine certain technical paramefers fo be im¬ 
plemented in our Monte Carlo code: Ihickness of fhe subsfrafe, defecfion zone deplh and palfern. 
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Figure 11. Left-hand plot: optical image and scanning electron microscope image (inset) of a cross-section 
of a S12571-100C MPPC showing the boundary between two pixels. Right-hand plot: optical image of a 
cross-section of a S13360-1350CS MPPC showing a pixel surrounded by optical trenches. The p-n junction 
and the thin p-H- electric contacts of the pixels can be distinguished. The trenches of the S13360-1350CS 
MPPC extend down to the bottom edge of the avalanche layer to prevent direct optical crosstalk. 


width and depth of the integrated wells, etc. Examples of these cross-sections for the S12571-100C 
and S13360-1350CS MPPCs are in figure 

All the MPPC series have a substrate with a thickness of around 630 /im and a similar layer 
structure of pixels. The p-type layer has an estimated thickness of 2 /rm. Between it and the 
substrate, there is a thin avalanche layer with thickness Tav ~ 1 /im, which we guessed to be lightly 
doped n type. In each pixel, a p-i-i- contact (~ 0.5 /im thickness) is on the top of the p-type layer. 
The coating and passivating surface layers have an estimated total thickness of 1.5 /f m. The S13360 
series have in addition trenches between pixels. As can be seen in figure o, fhe frenches exfend 
from fhe fop edge of fhe p layer down fo fhe boffom edge of fhe avalanche layer. The parficular 
pixel geomefry of each MPPC was implemenfed in fhe Monfe Carlo code. 

5.2 Algorithm 

The simulafion assumes fhaf fhe primary avalanche lakes place in a cenfral pixel of fhe array. A 
predefined number of secondary phofons for each discrete wavelenglh are generated in proportion 
lo fhe above-menlioned relalive speclrum. The inilial phofon posilion is uniformly disfributed in 
fhe volume of fhe pixel avalanche region, which is defined by a recfangular box of dimensions 
Lav X Lav X Lay. We assumcd Iwo possible values of fhe lateral size: Lav = 2.pix and Lav = Lpix • Vf, 
where F is fhe fill factor. The emission direcfion of fhe phofon is isolropically disfributed and 
ifs fofal dislance fraveled is generafed wilh a negative exponenfial disfribulion wilh paramefer j, 
where A is fhe absorption lengfh of silicon af fhe given wavelenglh. All fhe layers of fhe device, 
including fhe coating layer, are assumed fo have fhe same absorption lengfh. 

If a phofon is absorbed in fhe acfive region of a pixel differenl from fhe primary one, if Iriggers 
a prompf-crossfalk avalanche Ihere. The acfive region is defined by a recfangular box wilh fhe same 
lateral size as fhe avalanche region Lav and wilh a fhickness comprising bolh fhe avalanche layer and 
fhe p layer, buf nol fhe p-i-i- confacf. A 100% trigger probability is assumed, since the simulation 
is only intended to obtain relative results on the prompt crosstalk probability. Specular reflections 
on the top and bottom surfaces of the device are also simulated. For simplicity, both surfaces are 
assumed to be perfectly smooth and to have 100% reflectivity, regardless the photon incident angle 
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and wavelength. The reflected photons are tagged so that their contribution to crosstalk can be 
evaluated separately. 

For the SI3360 series, if the photon trajectory intersects a trench, the photon is dumped and 
the simulation continues to the next one. This means that direct crosstalk is not allowed taking into 
account the depth of the trenches (right-hand plot of figure [^. On the other hand, we considered 
no obstacle for transmission of photons through the coating layer above a trench. Therefore, the 
trenches do not eliminate completely the crosstalk due to photons reflected on the top-surface of 
the device, as will be shown below. 

If a photon is absorbed in the substrate, it produces an electron-hole pair. The diffusion of 
the minority caiTier, i.e., the hole, is simulated as a random walk process with an step of 0.1 /xm, 
but instead of forcing it to move in a discrete grid, the direction of movement in each step is 
randomly generated with a uniform distribution. The elapsed time is calculated as the number of 
steps times , where Dh is the diffusion constant of holes in the silicon substrate. The hole 

keeps diffusing until it reaches a border of the substrate or the total elapsed time is longer than 
300 ns, but recombination is ignored (i.e., Tbuik —^ °°)- If the hole reaches the active region of a 
pixel, it produces instantaneously a secondary avalanche and the event is classified as afterpulsing 
or crosstalk, depending on whether the pixel is the primary one or not. Crosstalk events with delay 
time shorter than 3 ns are classified as prompt crosstalk. The trigger efficiency of holes is also 
assumed to be 100% and the recovery of the primary pixel is ignored. Generation and transport of 
carriers in other parts of the device are not included either. 


5.3 Components of correlated noise 

In figure [T^ we show the wavelength distributions of the simulated photons causing the different 
components of afterpulsing and crosstalk for the S12571-050C MPPC, which has a pixel pitch of 
50 /im. The results were obtained assuming Dh = 0.4jUm^/ns and Lav = 39.4/im (i.e., F = 62%). 
The wavelength distributions for all the crosstalk components are within the range 600 - 1200 nm. 
This is determined by the absorption length of silicon, which varies from a few /rm at a wavelength 
of 600 nm to around 100 pm at 1000 nm and thereafter rises steeply 10. Photons with wavelength 
larger than 1200 nm are very unlikely to be absorbed within the device, and those with wavelength 
shorter than 600 nm do not scape the primary pixel or are absorbed very close to it, having a 
negligible probability to cause crosstalk. On the other hand, if assuming Lav = 50pm, i.e., no 
“dead zone” between pixels, the number of crosstalk events increases approximately by a factor 
of 4 and its wavelength distribution extends to shorter wavelengths. A similar but more moderate 
effect is observed when decreasing the pixel size. Apart from that, the wavelength range of photons 
contributing to crosstalk barely changes when varying other input parameters of the simulation. For 
instance, we obtained similar wavelength distributions for the crosstalk components when assuming 
that the photon emission spectrum is that of a black body of temperature 2000 K or 4500 K, as was 
also observed in [HBl . 

The prompt crosstalk due to top-side reflected photons is almost as high as the direct crosstalk, 
whereas the contribution of bottom-side reflected photons is negligible. Obviously, photon reflec¬ 
tions are highly overestimated in the simulation, as we assumed 100% reflectivity. However, it is 
worth considering the top-side reflections in the design of SiPMs, e.g., the metal circuitry and the 
trenches. In fact, for the simplified geometry implemented in the simulation for the S13360 series. 


- 18 - 




0.1 


—Photon emission spectrum 

- -AP (diffusion) 

- -CT (diffusion > 3 ns) 

— CT (diffusion < 3 ns) 

— CT (direct) 

=CT (top-side reflection) 
^■CT (bottom-side reflection) 



400 600 800 1000 1200 1400 1600 1800 

Wavelength (nm) 


Figure 12. Photon emission spectrum [@] used in the simulation and the wavelength distributions of the 
different components of crosstalk and afterpulsing obtained for the S12571-050C MPPC, assuming Lav = 
39.4/rm. Results are for a statistics of 10^ simulated photons. 

we observed that the trenches reduce the crosstalk component due to top-side reflected photons 
approximately by a factor of 4, but this would be completely eliminated if the baiTiers extend up to 
the top surface of the device. 

The simulation also overestimates the hole-induced afterpulsing and crosstalk, because it as¬ 
sumes an infinite bulk lifetime and 100% trigger probability for holes. This also distorts the shape 
of the wavelength distributions for these components, because the photon wavelength is correlated 
with both the depth at which the photon is absorbed and the time needed by the hole to reach a 
pixel. Most of the simulated holes causing afterpulsing and crosstalk are generated at a depth less 
than 10 /I m with respect to the bottom edge of the avalanche layer, although their depth distribution 
extends to several tens of /im. This distribution becomes flatter if assuming a diffusion constant of 
1 /fm^/ns instead of 0.4/rm^/ns. The holes contributing to prompt crosstalk (i.e., delay time shorter 
than 3 ns) are those generated in a thin layer of 3 /im beneath the active region of pixels different to 
the primary one. We found that the trenches of the S13360 series reduce the hole-induced crosstalk 
only by a few percent, with the reduction in the prompt component being smaller than 10%. Deeper 
trenches may reduce the crosstalk to a much greater extent. 

Unlike crosstalk, the wavelength distribution for afterpulsing in figure |T^ decreases with de¬ 
creasing wavelength approximately as the same rate as the photon emission spectrum does, because 
short-wavelength photons produce holes very close to the primary pixel. However, this distribution 
was obtained ignoring the recovery of the primary pixel. Including this would preferentially reduce 
the contribution of short-wavelength photons. On the other hand, the pixel size and the Lav param¬ 
eter were found to have a minor influence on the total number of afteipulses per primary avalanche. 
The impact of trenches is also very small. 

5.4 Spatial distribution of crosstalk 

We obtained for each crosstalk component the 2D histogram of fired pixels for a primary avalanche 
triggered in the central pixel of the array. In figure |TJ, the spatial distribution of direct crosstalk 
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Figure 13. Distribution of fired pixels via direct crosstalk for the S12571-100C (left) and S12571-010C 
(right) MPPCs. The simulation assumes that the primary avalanche takes place in the central pixel and does 
not include crosstalk cascades. The highlighted central square in each plot encloses 25 pixels. Results were 
obtained for a statistics of 10^ simulated secondary photons. 


obtained for the 100 and 10 /fm pixel pitch MPPCs of the S12571 series are shown, where Lav 
was assumed to be 88.3 and 5.74/im, respectively. Most of the crosstalk avalanches are triggered 
in the 24 pixels surrounding the primary pixel (i.e., the highlighted central square in each plot). 
In particular, the fraction of avalanches in this area is 92% (82%) for the S12571-100C (-010C) 
MPPC, with 84% (67%) of total being in the 8 nearest pixels and 71% (50%) in the 4 nearest 
ones. Assuming Lav = Cpix leads to the fractions 97% (89%), 93% (78%) and 85% (62%), in the 
same order. As observed with experimental data in section the larger the pixel size, the lower 
the number of neighboring pixels subject to direct crosstalk. Similar results were obtained for the 
prompt crosstalk component due to top-side reflected photons. On the other hand, for the SI3360 
MPPCs, the trenches eliminate the direct crosstalk and only allow photons with small incident 
angle to be reflected on the top surface. As a consequence, the totality of avalanches due to top¬ 
side reflected photons take place in the 8 nearest pixels. 

The spatial distribution for hole-induced crosstalk is more spread out. For instance, for the 
S12571-010C MPPC, the fraction of avalanches with delay time shorter than 300 ns in the 24 
nearest neighbors of the primary pixel is only 63%, and the fraction increases to 69% if the delay 
time is limited to 3 ns, i.e., the prompt component. 

Using these simulation results, we calculated the prompt crosstalk multiplicity distribution 
to be compared directly with the experimental data described in section This was caiTied 
out by means of a Monte Carlo algorithm that generates a large number of random cascades of 
prompt crosstalk, each one initiated by firing a random pixel of the MPPC. During a cascade, the 
probability that a fired pixel induces prompt crosstalk in other pixel (not fired yet) is proportional to 
the counts registered in the corresponding relative position in the above 2D histogram obtained by 
simulation. The proportionality constant is determined so that the sum of the crosstalk probabilities 
of all the pixels gives e. Actually, this sum depends on the position of the primary pixel in the 
array. Thus, we calculated it for all the possible primary pixel positions and then made an average. 
The S12571-050C MPPC was chosen for this analysis, since we measured an amplitude spectrum 
with many peaks and good resolution for it (figure Only the 2D histogram of simulated direct 
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Figure 14. Comparison of the prompt crosstalk multiplicity distribution calculated using simulation results 
of direct crosstalk for the S12571-050C MPPC with the experimental distribution obtained from the pulse 
amplitude spectrum of hgure The distributions calculated for three simple geometrical assumptions (i.e., 
the prompt crosstalk can only take place in the 4, 8 or 24 nearest neighbors of a fired pixel, with the same 
probability for all these neighbors) are also included in the comparison. 


crosstalk events was used in the calculation, because the top-side reflected photons lead to a very 
similar spatial distribution, as explained above, and the contribution of holes to prompt crosstalk 
is supposed to be small for this MPPC. For simplicity, this 2D histogram was averaged over all 
the pixels outside the square enclosing the primary pixel and its 24 nearest neighbors, that is, we 
assumed that a fired pixel has the same probability to induce prompt crosstalk in any pixel that is 
not within this area. This approximation can be done because only 10% of the simulated prompt 
crosstalk avalanches take place outside this area of pixels for the S12571-050C MPPC. We assumed 
Lav = 39.4/im (i.e., F = 62%) in the simulation. 

The relative difference of the calculated prompt crosstalk multiplicity with respect to the exper¬ 
imental one for the S12571-050C MPPC is represented by blue circles in figure [T^ A remarkable 
agreement is observed, although the difference for some multiplicities is larger than the error bars, 
which include both the experimental uncertainties and the statistical uncertainties of the numerical 
calculation. We also show the comparison for three multiplicity distributions calculated in the same 
way, but instead of using the simulation results, assuming that the prompt crosstalk can only take 
place respectively in the 4, 8 or 24 nearest neighbors of a fired pixel and that all these neighbors 
have the same crosstalk probability. Unlike the analytical formulae used in section ^3L these nu¬ 
merical calculations allow the determination of the distribution for arbitrarily high multiplicities 
and take into account that border pixels have a smaller number of neighbors subject to prompt 
crosstalk. 

5.5 Time distribution of afterpulsing and delayed crosstalk 

Our simulation also provides the diffusion time distribution of simulated holes that cause after¬ 
pulsing and crosstalk. Results for the S12571-050C MPPC, assuming Lav = 39.4jUm and = 
0.4/im^/ns, are shown in figure |T^. A power law fit to data between 3 and 100 ns is shown for both 
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Figure 15. Diffusion time distribution of simulated holes that cause afterpulsing (red squares) and crosstalk 
(blue squares) for the S12571-050C MPPC. Solid lines represent power law fits. Filled circles represent 
the time distribution after including the effects of the primary pixel recovery and hole recombination in the 
substrate, assuming Tree = 20 ns and Tbuik = 10 ns. 


components, with the fitted exponent being a = —0.96 for afterpulsing and b = —0.50 for delayed 
crosstalk. Similar a and b values are obtained for other simulation input data, where a ranges from 
-0.83 to -1.25 and b ranges from -0.24 to -0.56. A slight decrease in \a\ as increasing pixel size or 
decreasing diffusion constant is observed. The dependence of b on these parameters is not clear, 
but trenches make the time distribution for crosstalk flatter, with b > —0.30, because they intercept 
some laterally emitted photons that would generate holes close to the avalanche layer (i.e., short 
diffusion time). In general, power law functions with exponents a = —1.0 and b = —0.5 fit well 
to the simulation data and also allow us to describe the experimental delay time distribution of 
secondary pulses, as described in section 

As mentioned above, neither the recovery of the primary pixel or the recombination of holes 
in the substrate are simulated. In figure |T^, we also show the time distributions of afteipulsing and 
crosstalk after including these effects as described in section EH where Tree = 20 ns and Tbuik = 
10 ns are assumed. 

We calculated the relative contributions of afterpulsing and delayed crosstalk as a function 
of the pixel size. To do that, we refitted the simulation data assuming a = —1.0 and b = —0.5. 
The following ^ ratios were obtained for the MPPCs of the S12571 series: 0.14 for Lp^ = 
10/fm, 0.030 for Lpix = 25/fm, 0.013 for Lpix = 50/fm and 0.006 for Lpix = 100/im, which are in 


general agreement with the ratios obtained experimentally in section 1.2 


The decrease in ^ with 

Cap 


increasing pixel size is due to the fact that short-wavelength photons contribute differently to the 
afteipulsing and crosstalk components, as explained in section p3| . The results are only slightly 
dependent on the input parameters of the simulation, e.g.. Lav and Db. 


6. Conclusions 

The correlated noise in several Hamamatsu MPPCs was analyzed in depth. We considered the noise 
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contribution of the secondary photons emitted by avalanches as well as of the photon-generated 
minority carriers diffusing in the silicon substrate of the device. Our experimental method was 
based on a waveform analysis that is able to identify pulses with time differences as low as 3 ns. 
This determined our empirical criterium to differentiate between prompt and delayed crosstalk. 
In addition, we used the fact that the afterpulse amplitude is lowered during the primary pixel 
recovering to characterize the afterpulsing separately from the delayed crosstalk. 

The prompt crosstalk was characterized experimentally through the pulse amplitude spectrum 
at dark conditions. In particular, the gain dependence of the prompt crosstalk probability was mea¬ 
sured and parameterized. The new S13360 series of MPPCs, which incorporate isolating trenches 
between pixels, were found to have a prompt crosstalk probability more than a factor of 30 lower 
than the former MPPC series. However, it is still non-negligible, likely due to both secondary 
photons reflected on the top surface of the MPPC and minority carriers diffusing in the substrate. 
Using the information provided by the prompt crosstalk multiplicity distribution and analytical 
models proposed in a previous work [0, we deduced that prompt crosstalk mostly takes place in a 
small area of pixels (~ 8) around the primary one. 

A novel procedure to characterize the afterpulsing and delayed crosstalk was presented in this 
work. It uses both the amplitude and delay time distributions of secondary pulses. The primary 
pixel recovery, the pulse detection threshold and the prompt crosstalk in secondary avalanches were 
carefully taken into account. In addition, the time distributions of both afterpulses and delayed- 
crosstalk pulses due to diffused minority carriers were modeled using results from a Monte Carlo 
simulation. A significant component of delayed crosstalk was observed in all the MPPCs that 
were characterized. Also, our results indicate that the dominant contribution of afterpulsing is 
due to diffused minority carriers, not to the release of avalanche carriers trapped in deep levels 
of the depleted layer, as commonly believed. In fact, we attributed the drastic afterpulse reduction 
observed in the new MPPC series to an increase in the doping density of the n-type silicon substrate 
from 3-6-10*^ cm^^ to 2 • 10^^ cm^^, involving a notably decrease of the minority carrier lifetime 
from a few hundreds of ns to around 10 ns. The recovery time was also shown to have an important 
role in the afterpulse reduction. 

We developed a Monte Carlo simulation of the physical processes causing afterpulsing and 
crosstalk. It includes the emission and transport of secondary photons as well as the diffusion of 
photon-generated minority carriers (holes) in the substrate. We measured the geometrical param¬ 
eters of the various MPPCs to be implemented in the code. The simulation allowed us to analyze 
in detail the different components of afteipulsing and crosstalk, e.g., the wavelength distribution 
of secondary photons related to each noise component, the spatial distribution of crosstalk events, 
and the diffusion time distribution of minority carriers causing afterpulsing and delayed crosstalk. 
The simulation results are in general agreement with experimental data and support our hypotheses 
for their interpretation. We also concluded that, although the trench structure used in the S13360 
MPPCs proved to reduce drastically the prompt crosstalk, longer trenches extending from the top 
surface of the device down to several /rm into the substrate would reduce both prompt and delayed 
crosstalk to a much greater extent. 
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